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'nI" ■ 1 Introduction 
cn 

^' 

I group-invariant (called also similarity) solutions of partial differential equations 

■ (PDEs) [3]. It was called by the authors "non-classical" to emphasize the differ- 
ence between it and the "classical" Lie reduction method described, e.g., in [16,17]. 

\ A precise and rigorous definition of nonclassical invariance was firstly formulated 

^ ■ in [11] as "a generalization of the Lie definition of invariance" (see also [26]). Later 



As early as in 1969 Bluman and Cole introduced a new method for finding 



operators satisfying the nonclassical invariance criterion were called, by different 
authors, nonclassical symmetries, conditional symmetries and Q-conditional sym- 
metries [8,10,15]. Until now all names are in use. Following [19] we call nonclas- 
sical symmetries reduction operators. The necessary definitions, including ones of 
equivalence of reduction operators, and relevant statements on this subject are 
collected in [22]. 

The problem of finding reduction operators for a PDE is more complicated 
than the similar problem on Lie symmetries because the first problem is reduced 
to the integration of an overdetermined system of nonlinear PDEs, whereas in 
the case of Lie symmetries one deals with a more overdetermined system of linear 
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PDEs. The complexity increases in times in the case of classification problem of 
reduction operators for a class of PDEs having nonconstant arbitrary elements. 

Often the usage of equivalence and gauging transformations can essentially 
simplify the group classification problem. Moreover, their implementation can 
appear to be a crucial point in solving the problem. This observation is justified 
by a number of examples [12,21,22]. The above transformations are of major 
importance for studying reduction operators since under their classification one 
needs to surmount much more essential obstacles then those arising under the 
classification of Lie symmetries. 

In [22] simultaneous usage of equivalence transformations and mappings be- 
tween classes allowed us to carry out group classification of the class of variable 
coefficient semilinear reaction-diffusion equations with power nonlinearity 

f{x)ut = {g{x)u^)^ + h{x)u"', (1) 

where / = f{x), g = g{x) and h = h{x) are arbitrary smooth functions of the 
variable x, f{x)g{x)h{x) 7^ 0, m is an arbitrary constant (m / 0, 1). 

In the same paper an algorithm for finding reduction operators of class (1) 
involving mapping between classes was proposed. Here using this algorithm we 
investigate reduction operators of the equations from class (1) with m / 2. The 
case m = 2 will not be systematically considered since it is singular from the Lie 
symmetry point of view and needs an additional mapping between classes (see [22] 
for more details). Nevertheless, all the reduction operators constructed for the 
general case m 7^ 0, 1, 2 are also fit for the values m = 0, 1, 2. 

The structure of this paper is as follows. For convenience of readers sections 
2-4 contain a short review of results obtained in [22] and used here. Namely, in 
section 2 all necessary information concerning equivalence transformations and 
mapping of class (1) to the so-called "imaged" class is collected. Results on group 
classification and additional equivalence transformations of the imaged class are 
also presented. Section 3 describes the algorithm for finding reduction operators of 
class (1) using mapping between classes. In section 4 known reduction operators 
of constant-coefficient equations from the imaged class are considered. Their 
preimages are obtained. The results of sections 5 and 6 are completely original 
and concern the investigation of reduction operators for equations from the imaged 
class which have at least one nonconstant arbitrary clement. It appears that 
application of the reduction method to equations from the imaged class with 
m = 3 leads in some cases to necessity of solving first-order nonlinear ODEs of 
a special form related to Jacobian elliptic functions. The table with solutions of 
the ordinary differential equations (ODEs) of this kind is placed in Appendix. 

2 Lie symmetries and equivalence transformations 

To produce group classification of class (1), it is necessary to gauge arbitrary 
elements of this class with equivalence transformations and subsequent mapping 
of it to a simpler class [22] . 
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Theorem 1. The generalized extended equivalence group of class (1) is formed 
by the transformations 

t = dit + d2, x = (p{x), u = i/j{x)u, 

f = 7^ J J 9= — To" ff) = 1 — zt"; m = m, 

where ip is an arbitrary smooth function of x, (fx ^ Q and ip is determined by the 
formula ij{x) = {S3 J + (^4) . Sj (j = 0,1,2,3,4) are arbitrary constants, 

SoSi{S3^ + S^)^0. 

The usual equivalence group G~ of class (1) is the subgroup of the generalized 
extended equivalence group G"", which is singled out with the condition S3 = 0. 

The presence of the arbitrary function (p{x) in the equivalence transformations 
allows us to simplify the group classification problem of class (1) via reducing 
the number of arbitrary elements and making its more convenient for mapping to 
another class. 

Thus, the transformation from the equivalence group C 

t' = sign{f{x)gix))t, x' = j \f{x)g{x)-^\^ dx, u' = u (2) 



maps class (1) onto its subclass f'{x')u'ti = {f {x')u' x')x' + h'{x')u'"^ with the 
new arbitrary elements m' = m, f'{x') = g'{x') = sign(5(a;)) |/(.t)(/(.t)| 2 and 

h'{x') = \g{x)f{x)^^\^ h{x). Without loss of generality, we can restrict ourselves 
to study the class 

fix)ut = {fix)u,)x + h{x)u"', (3) 

since all results on symmetries and exact solutions for this class can be extended 
to class (1) with transformation (2). 

It is easy to deduce the generalized extended equivalence group for class (3) 
from theorem 1 by setting f = g and f = g. See theorem 4 in [22]. 

The next step is to make the change of the dependent variable 



vit,x) = ./\f{xj\u{t,x) (4) 

in class (3). As a result, we obtain the class of related equations of the form 

vt = v^^ + H{x)v'^ + F{x)v, (5) 

where the new arbitrary elements F and H are connected with the old ones via 
the formulas 



F(a;) = -i^^b, i/(x) = ^^^^iSM. (6) 

Since class (5) is an image of class (3) with respect to the family of transforma- 
tions (4) parameterized by the arbitrary element /, we call them the imaged class 
and the initial class, respectively. 
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Theorem 2. The generalized extended equivalence group Gpjj of class (5) coin- 
cides with the usual equivalence group Gpjj of the same class and is formed by 

the transformations 



i=SiH + 62, 



X 



F = 



Si'' 
h 



H = 



5lX + ^3, 

H 



64V, 



X 2 X m—1 ' 



m = m. 



where Sj, j = 1, ... ,4, are arbitrary constants, S164 7^ 0. 

The following important proposition is proved in [22]. 

Proposition 1. The group classification in class (1) with respect to its generalized 
extended equivalence group G~ is equivalent to the group classification in class (5) 
with respect to the usual equivalence group Gpjj- of this class. A classification list 
for class (1) can be obtained from a classification list for class (5) by m,eans of 
taking a single preimage for each element of the latter list with respect to the 
resulting mapping from class (1) onto class (5). 

All possible G^ji^-inequivalent values of the parameter-functions F and H for 
which equations (5) admit extension of Lie symmetry are listed in table 1 together 
with bases of the corresponding maximal Lie invariance algebras. 

Table 1. The group classification of the class vt = Vxx + H{x)v"^ + F{x)v. 

m 7^ 0.1,2; H(x) 0. 



N 


H{x) 


F{x) 


Basis of A"^^ 





V 


V 


dt 


1 




ai 


dt, dx + avdy 


2 




-a^ 


dt, dx + avdy, 2tdt + (x — 2at)dx+ 
(a(x - 2at) + ^) vd^ 


3 


5x^ 


a2X~'^ 


dt, 2tdt + xdx + ^ vd^ 


4 


6x^6^^^ 


-p'^x^ + a2X-2 + 7 


dt, e^^^dt + 2l3xe^^*dx- 
2(3e^f^'{px^ ^^l)vd. 


5 






dt, e2^* [dx - Pxvd^] 


6 






dt, e^^* [dx - fixvd^] , 

6^/3* ^dt + 2(5xdx - 2/3 - ^) vd^\ 



Here a, /?, 7, 5, k,p, q, ai , 02, 03 are constants satisfying the conditions: a = 
f3= ^,1 = S = ±lmodGpjj; p ^ 0, ^ -a\ ^ + ^ 0, 

q'+al^O, a3/fe^. 
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The results on group classification of class (3) can be found in table 3 of [22]. 

Additional equivalence transformations between G^^-inequivalent cases of Lie 
symmetry extension are also constructed. The independent pairs of point-equivalent 
cases from table 1 and the corresponding transformations are exhausted by the 
following: 



6 1-^ 2|g=o '■ the previous transformation with k = 0. 

The whole set of form-preserving [13] (also called admissible [18]) transformations 
of the imaged class for the case m 7^ 0, 1, 2 is described in [22]. 

3 Construction of reduction operators using 
mappings between classes 

Here we adduce the algorithm of application of equivalence transformations, gaug- 
ing of arbitrary elements and mappings between classes of equations to classifica- 
tion of reduction operators. 

1. Similarly to the group classification, at first we gauge class (1) to subclass (3) 
constrained by the condition f = g. Then class (3) is mapped to the imaged 
class (5) by transformation (4). 

2. Since nonclassical symmetries of constant coefficient equations from the im- 
aged classes are well investigated (see below for more details), they should be 
excluded from the consideration. It also concerns variable coefficient equations 
from class (5) which are point-equivalent to constant coefficient ones, namely 
equations associated with cases l|g^O) 2|q^o and 6 of table 1 and equations re- 
duced to them by transformations from the corresponding equivalence groups. As 
a result, only equations from class (5) which are inequivalent with respect to all 
point transformations to constant coefficient ones should be studied. 

3. Reduction operators should be classified up to the equivalence relations 
generated by the equivalence group or even by the whole set of admissible trans- 
formations. Only the nonsingular case r ^ (reduced to the case r = 1) should 
be considered. Operators equivalent to Lie symmetry ones should be neglected. 

4. Preimages of the obtained nonclassical symmetries and of equations admit- 
ting them should be found using backward gauging transformations and mappings 
induced by these transformations on the sets of operators. 

Reduction operators of equations from class (3) are easily found from reduction 
operators of corresponding equations from (5) using the formula 



1 1 



q=0, ai=ai+a^ i 



2^2 



q=o ■ t = t, X = X + 2at, V = e ; 





(8) 
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Here r, ^ and rj are coefficients of dt, dx and respectively, in the reduction 
operators of equations from class (5). The substitution v = \/\f\u is assumed. 

There exist two ways to use mappings between classes of equations in the in- 
vestigation of nonclassical symmetries. Suppose that nonclassical symmetries of 
equations from the imaged class are known. The first way is to take the preimages 
of both the constructed operators and the equations possessing them. Then we 
can reduce the preimaged equations with respect to the corresponding preimaged 
operators to find non-Lie solutions of equations from the initial class. The above 
way seems to be non-optimal since the ultimate goal of the investigation of non- 
classical symmetries is the construction of exact solutions. This observation is 
confirmed by the fact that the equations from the imaged class and the associated 
nonclassical symmetry operators have, as a rule, a simpler form and therefore, are 
more suitable than their preimages. Reduced equations associated with equations 
from the imaged class are also simpler to be integrated. Moreover, it happens 
that preimages of uniformly parameterized similar equations do not have similar 
forms and belong to different parameterized families. As a result, making re- 
ductions in the initial class, we have to deal with a number of different ansatze 
and reduced equations although this is equivalent to the consideration of a single 
ansatz and the corresponding reduced equation within the imaged classes. This 
is why the second way based on the implementation of reductions in the imaged 
classes and preimaging of the obtained exact solutions instead of preimaging the 
corresponding reduction operators is preferable. 



4 The case of constant F and H 

Constant coefficient equations from the imaged class belong to the wider class of 
a quasilinear heat equations with a source of the general form vt = Vxx + Qiv)- 
Lie and nonclassical symmetries of these equations were investigated in [5,6] and 
[2,4,9,20], respectively. Their non-Lie exact solutions were constructed by the 
reduction method in [2,4], see also their collection in [22]. The nonlinear equation 
= Vxx + q{v) possesses pure nonclassical symmetry operators with nonvanishing 
coefficients of dt if and only if g is a cubic polynomial in v. Thus, in the case 
q = Sv"^ + ev, where S 0, such operators arc exhausted, up to the equivalence 
with respect to the corresponding Lie symmetry groups, by the following: 

S<0: dt±^V^vdx + l{Sv^ + ev)dy, 
s = 0: dt-^dx-^vd^, 
£ < 0: dt + 3/xtan(/xx)5x — Sfi^ sec^{iJ,x)vdv, 
£ > 0: dt — 3iJ,ieiXi.h.{iJ,x)dx + 3/x^sech^(/Ltx)v5„, 
dt — 3/xcoth(/xx)5a; — 3/x^cosech^(/Ltx)ti5^, 



where /x = -^|£|/2. Note that the last operator was missed in [2,4]. 
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Finding the preimages of equations with such values of q with respect to trans- 
formation (4) and the preimages of the corresponding reduction operators ac- 
cording to formula (8), we obtain the cases presented in tabic 2. In this table 
cl + C2 0, 1^ > 0, n = I'/V^; £ = 0, £ = i''^ > and e = -i^'^ < in cases 1, 2 
and 3, respectively. 



Table 2. Nonclassical symmetries of equations of the form 
f{x)ut = {f{x)u,), + 6f{xfu\ fix) = C(x)2. 



N 


C(^) 


Reduction operators 


1 


CiX + C2 


Ot Ox TT'^^u 

x x^C, 








2 


ci sm{i'x)+ 

C2 COsiux) 


dt — 3/Utanh(/ia;)9x + 3/i tanh(/ia;) + /usech^(/ix)^ u9„, 

(9/ — 3// cothf/crlr),. + 3// ^— cotlif/rr) — // (•()S('cli^(//:r)^ udu 


3 


c\ sinh(i^a;)+ 


dt ± y-25Cudx + T V-2(5C^M - y^)udu, 




C2 cosh(i/a;) 


dt + 3^ta.n{ij,x)dx — 3/i tan(/xa;) + /isec2(/ia;)^ udu 



5 Reduction operators for general values of m 

In this section we look for G^ji^-inequivalent reduction operators of the imaged 
class (5). Here reduction operators have the general form Q = rdt + + rjdy, 
where r, ^ and rj are functions of t, x and v, and (r, ^) ^ (0,0). Since (5) is an 
evolution equation, there are two principally different cases of finding Q: t ^ 
and T = [10, 14,25]. The singular case r = was exhaustively investigated for 
general evolution equation in [14,25]. 

Consider the case t ^ 0. We can assume r = 1 up to the usual equivalence of 
reduction operators. Then the determining equations for the coefficients ^ and rf 
have the form 

Vt - Vxx + 2Cx?? = 

C {Exv"" + Fxv) + {2ix - Vv) {Hv"" + Fv)+r]{F + Hv^^-^m) , 

3^^ (Hv"" + Fv) + 2^^^ + it + '^Vvx - ^xx - 24^r/ = 0. 
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Integration of first two equations of system (9) gives us the following expressions 
for ^ and rj 

^ = av + b, 

1 (10) 

77 = —-a^v^ + {ux — ab)v^ + cv + d, 

where a = a{t, x), b = b{t, x), c = c{t, x) and d = d{t, x). 

Substituting ^ and r] from (10) into the third and forth equations of (9), we 
obtain the classifying equations which include both the residuary uncertainties in 
coefficients of the operator and the arbitrary elements of the class under consid- 
eration. 

Since the functions a, b, c, d, F and H do not depend on the variable v, the 
classifying equations should be split with respect to different powers of v. 

Two principally different cases a = and a 7^ should be considered separately. 

If a = then for any m 7^ 0, 1, 2 the splitting results in the system of five 
equations 

mHd = 0, dt-dxx + '^bxd - Fd = 0, 
bt - bxx + ^bbx + 2cx = 0, 
bHx + (c(m - 1) + 2bx) H = 0, 
bFx + 2bxF + Cxx — ct — 2bxC = 0. 

Since mH 7^ then d = and the second equation of (11) becomes identity. 

Finding the general solution of the other three equations from (11) appears to 
be a very difficult problem. But it is easy to construct certain particular solutions 
setting, e.g., bt = 0. This supposition implies that q = 0. Then the integration 
of (11) gives the expressions of c, F and H via the function b{x) 

1 2 1 
c = --6 +-^bx + ki, 



where ki, k2 and ^3 are arbitrary constants, ^3 7^ 0. 

Theorem 3. The equations from class (5) with the arbitrary elements given by 
formulas (12) and (13) admit reduction operators of the form 



where b = b{x) is an arbitrary smooth function and ki is an arbitrary constant. 




(12) 



(13) 




(14) 
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Note 1. Theorem 3 is true for any m G M, including m G {0, 1, 2}. 

We present illustrative examples, by considering various forms of the function b. 

Example 1. We take b = x"^ and substitute it in formulas (12)- (14) to find that 
the equations 



= + _^gi(m-l)(.3+6fci.-i)^m + (Aa^i + 2x + ki)v, (15) 



admit the reduction operator 

Q = dt + x^dx + + X + ki^vdy. 

The corresponding ansatz v = xe^'^*~&^^z{u), where u = t + x~^, gives the 
reduced ODE 

2u;u, + A;3e'=i("*-^)'"z™ + /c2^ = 0. (16) 
For fci = the general solution of (16) is written in the implicit form 

r 1 f 2fc3 ^m+l ^ ^ _i 

{Z-k2Z^ + C^)-^dz = ±uj + C2, Z={ ' '^^ (17^ 

J I— 2A;3ln2;, m = — 1. 

If /c2 = and m ^ —1, we are able to integrate (17). Setting Ci = 0, we obtain 
a partial solution of the reduced equation in an explicit form: 

2 

1— m 



where C is an arbitrary constant. Note that the constant C can be canceled via 
translations of uj induced by translations of t in the initial variables. 
In the case A;i = A;2 = and m / — 1 we construct the exact solution 



V = xe 



2 

1— m 



of the corresponding equation (15). Preimages of them with respect to transfor- 
mation (4) are the equation 



e~3^''x^'Ut = {e'^'^^ x^Ux)x + h x~'^e''^^^u"^ (19) 
and its exact solution 



u 



2 

1— m 
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Analogously, if /c2 = Ci = and m = —1 then integration of (17) gives 

erf (V — In zj = iby^^ w + C, where erf(y) = e~^^dt is the error func- 

tion and C is an arbitrary constant which can be canceled by translations of lo. 
Therefore, 



z = exp 



erf 



(20) 



where erf ^ is the inverse error function, represented by the series 

2k+l 



/ / — \ ■ 



, where cq = 1, 



Cfe 



m=0 



CmCk-l-m ^ li I 1 ^ 

(m + l)(2m + l) 1 ' ' 6' 90 ' 



The corresponding exact solution of equation (19) with m = — 1 is 

2"! 



u = exp 



Example 2. Consider b = x ^. In view of theorem 3 the equations from class (5) 
with the arbitrary elements 

F = ki+ k2X^ - 2x-2, H = fe3a;"^+^e^(i-™)'=i^' (21) 

admit the reduction operator 

Q = dt + x~^dx + (fci - x""^) vdy. 

The ansatz constructed with this operator is f = x~^e''^^z{Lo), where u = 
x^ — 2t, and the reduced equation reads 

4^^u, + A:3e^(i-'")'=i^z"* + k2Z = 0. 

If ki = k2 = 0, the reduced equation is integrated analogously to equation (16) 
and has the similar particular solution 



z = 



. m— 1 



2(m+l) 



UJ 



2 

1 — m 



exp 



erf-i ( ±^J^uj 



m 7^ — 1, 
-1. 



Substituting the obtained z to the ansatz, wc construct exact solutions of equa- 
tions from class (5) with the arbitrary elements (21) for the values ki = k2 = 0. 
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The preimaged equation x'^ut = {x^Ux)x + ^3 x^^'^'^^^'^ has the exact solution 



-3 / im-1 



X exp 



erf 



f ix' - 2t) 



m = —1. 



In the following two examples wc assume that fci = A;2 = in the formulas (12)— 
(14) since this supposition allows us to find preimages in class (3) with arbitrary- 
elements being elementary functions. 

Example 3. Let b = and ki = k2 = 0. The equations of the form 

vt = + A;3eM(i-™)^""+("^+3)^)t;"^ - 1 (e'^^ + 4e-^ + l)v, (22) 
admit the reduction operator 

Q = dt + e-^dx - ^ [e-^ + e-2-) vd^. 

The corresponding ansatz v = e^^^ ''~^^z{u), where u = e'^ — t, gives the 
reduced ODE 



z^^ + k^z^ = 0. 

It coincides with the equation (16) with k\ = k2 = 0, which has the particular 
solution (18) (rcsp. (20)) for m / — 1 (rcsp. m = —1). Substituting these solutions 
to the ansatz, we obtain exact solutions of equation (22). 

A prcimage of (22) with respect to transformation (4) is the equation 



e*" ""ut = (e^ '^Ux)x + k^e 
having the solution 



U 



2 

1— m 



cxpj- erf-i |^±y^(e^-t)^ |, m = -1. 

Example 4. Substituting 6 = sinx and /ci = ^2 = to formulas (12)-(14) and 
making the reduction procedure, we obtain the following results: The equation 

Vt = Vxx + kz{smx)~^^'^^''^^e^~^^°^^v'^ + ^ (cos^ x + 4cosa; + cosec^x)f , 
has the exact solution 



V = i 



6 2 '^°^^Vsinx (± 



Vsinx exp< 



m—l 



-a^(t-ln|tan| 



2 

1 — m 



COS x 

2 



erf-^ ±V^ (t-ln|tanf|) 



m + -1, 



, m = —1. 
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The corresponding equation from class (3) is 

gcos X gjj^ xut = {e^°^ ^ sin xux)x + ^3 cosec x e'^°^ ^u"^ 

whose exact solution is easy to be constructed from the above one using for- 
mula (4). 

We have shown the applicability of theorem 3 for construction of non-Lie exact 
solutions of equations from classes (5) and (3). Moreover, using these solutions 
one can find exact solutions for other equations from (5) and (3) with the help of 
equivalence transformations from the corresponding equivalence groups. 

Note 2. In the case m = 3 we are able to construct more exact solutions of 
equations from class (5) whose coefEcients are given by (12)-(13) with ki = 0, 
namely, for the equations 

= ... + fe35-e/--.3 +(^^-\b^ + b. + \ (^) ' - j (23) 

where b = b{x), / 0. 

According to theorem 3, equation (23) admits the reduction operator (14) (with 
ki = 0). An ansatz constructed with this operator has the form 



v = z{u)^\e-^^^'^'', where u = t- 




and reduces (23) to the second-order ODE 
Zuiui = —ksz^ — k2Z. 

It is interesting that the reduced ODE does not depend on the function b{x). 
Multiplying this equation by z^; and integrating once, we obtain the equation 

z^ = — -z - k2Z + Ci. 

Its general solution is expressed via Jacobian elliptic functions depending on values 
of the constants k2, k^ and Ci. See Appendix for more details. 

For example, if A;2 = 1 + /i^, k^ = —ijjp' and Ci = 1 (0 < < 1) we find two 
exact solutions of equation (23) 




where sn(a;, //), cd(a;,//) are Jacobian elliptic functions [24]. 

The second case to be considered is a 7^ 0. Then after substitution of ^ and r] 
from (10) to system (9) its last equation takes the form 
2 

-a^u^ -I- 2a{ab — 2a^)f;^ -|- (at -|- ^axx + 3aF - 2{ab)x — 2ac)v + ^^4^ 

bt + 2bxb - bxx - 2ad + 2c^ + ZaHv'^ = 0. 

It is easy to sec that a 7^ if and only if m = 3. The investigation of this case is 
the subject of the next section. 
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6 Specific reduction operators for the cubic 
nonlinearity 

Splitting equation (24) in the case m = 3 and a 7^ with respect to u, we obtain 
that the functions a, b c and d do not depend on the variable t and are expressed 
via the functions F and H in the following way 



2v^V 2i7 \h)^ 2\H 



(25) 



where £ = ±1. li H < Q the corresponding reduction operators have real coeffi- 
cients. 

Then splitting of the third equation of system (9) for m = 3 results in the 
system of two ordinary differential equations 

H^Hxxxx — 13 -f^;^ + 2FxH^Hx + 22 HH^Hxx — 4:FH^H^— 
4 H^H^^ - 6 H^H^H^^^ + 4 FH^H^^ - 6 F^^H"" = 0, 

IGF^a^^H^ + 16 H^H^H^^ + 3 H^H^H^^^ - 4 F^H^H^^- (26) 
6 H^HxxHxxx — 18 HH^Hxx — 8 FF^H^ + 2 FxH^H^— 
20 FH'^Hl - 12 FH'^Hxxx + 5 + 32 FH^HxHxx = 0. 

The following statement is true. 

Theorem 4. The equations from class (5) with m = 3 and the arbitrary elements 
satisfying system (26) admit reduction operators of the form 

Q = dt+ (^^V2eV^v + ^^dx + 



V2e / IHx {Hx\ 1 fHx 



2v^ \ " 2 H \H J ^ 2\H 
where e = ±1. 

Let us note that system (26) can be rewritten in the simpler form in terms of 
the functions F and b 

bxxx = 6Fxx-2 bxb"^ + bx^ + 2 bbxx - 2 - 4 Fbx, 

16 Fxxx = 4 bxFx + 2 bx'^b + 6 bxbxx + 2 b'^Fx + b%x+ (28) 

3 b'^bxx + 12 Fbxx + 8 FFx + 4 Fbbx ■ 
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System (26) consists of two nonlinear fourth- and third-order ODEs. Unfortu- 
nately we were not able to find its general solution. Nevertheless, we tested the 
six pairs of functions F and H appearing in table 1 in order to check whether they 
satisfy system (26). In the case of positive answer the corresponding reduction 
operator is easily constructed via formula (27). It appears that system (26) is 
satisfied by F and H from cases 1, 2 and 6 and by those from cases 3 and 4 for 
special values of the constants k and 02, namely, {k,a2) G |) , (— |, 

So, we can construct preimages of these equations using formulas (6). Below 
we list the pairs of the coefficients / and h for which the corresponding equations 
from class (3) with m = 3 admit nontrivial reduction operators. 

Hereafter bf + b2 7^ 0. The numbers of cases coincide with the numbers of 
the corresponding cases from table 1. (Case 5 does not appear below since the 
functions F and H from this case of table 1 do not satisfy system (26).) 

1. ai = 0: / = (6lX + 62)^ /i = (5e'^"(6ix + 62)^ q ^ 0; 
ai > ^ ai = 1 mod Gp^ : / = (&i sinx -|- 62 cos x)^, 

h = de'^^{bi sinx + 62 cosx)^; 
ai < ai = — 1 mod Gp^ : f = {bi sinhx + 62 coshx)^, 

^ = 5e'^^(6i sinhx + 62 coshx)"^, q ^ ±2. 

2. q = 0: /=(6lX + 62)^ /t = J(6ix + 62)^ 

q 7^ ^ q = —2 mod Gpjj : f = (bi sinhx + 62 cosh x)^, 

h = 6e~^^{bi sinhx -|- 62 coshx)'^. 



3. [k, 02) = (-3, 1) : / = x{bi sin(\/21n |x|) + 62 cos(y2 In |x|))^ 

h = fe-^(6i sin(\/21n|x|) + 62 cos(V2 In |x|))^ 
ik,a2) = (-f,^) : / = x(6i|x|3 +b2\x\-^f, 

h = d^/\x\{bi\x\^ + b2\x\~^)^. 



A.f = x-i (6iM«,^(px2) + b2W^,^{px'')f , 
h = Sx'^-'^eP^' (6iM«,^(px2) + b2W,,^{px^))^ 



where k = ^ = ^^V^, (A;,a2) G {(-3,|) , (-f,^)}. M«,^ and W«,^ 

are the Whittaker functions [24]. 



6. / = x~^ [biM_i 1 (px^) + b2W_i 1 (px^ 

\ 4 ' 4 4 ' 4 

h = dx-^eP""' (biM_i lipx^) + b2W_i 1 ipx^)\ 

\ 4'4 4'4 / 



Note that in the case p > the above Whittaker function is expressed via the 
error function: M_i i(px^) = 5 a/tt ^/px"^ ^ eii^s/px^) [1]. 
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Since the cases l|g^O) ^\q^o ^^'^ 6 are reduced to constant-coefficient ones we 
do not consider them. 

Example 5. Class (5) contains equations with cubic nonlinearity, which are not 
reduced to constant-coefficient ones by point transformations and admit reduction 

operators of the form (27). One of them is the equation with the coefficients F 
and H presented by case 3 of table 1 with k = —3, a2 = | and 6 = —1, namely, 

vt = Vxx - V + ^x-^v. (29) 
According to theorem 4 this equation admits two similar reduction operators 

(6 = ±1) 

Q± = dt + ^V2 (^ex-^v - V2x-^yx - 

'-^V2 (y2x~^v^ - Sex'iv'^ - V2x~\ + 4ex~^^dy. 



3 

4' 



They lead to the solutions differing only in their signs. Since equation (29) is 
invariant with respect to the transformation v 1— > —v, we consider in detail only 
the case e = 1. For all expressions to be correctly defined, we have to restrict 
ourself with values x > 0. (Another way is to replace x by \x\.) 
For convenient reduction we apply the hodograph transformation 

t = V, X = X, V = t 

which maps equation (29) and the reduction operator to the equation 

Vi^ Vxx + Vx"^ %- - 2 Vf V-x Vl^ + V-^ ^^'^t = ^ ^^'^^ 

and its reduction operator 

Q+ = -^\/2(^\/2x~^t^-3x-iP- V2x"2t-F4x-i)9j + 
^V2 (x-it - \f2x-^^dx + 5^5, 
respectively. An ansatz constructed with the operator (3+ has the form 

1^2^+V^ 1-2 u -2*-^^ 



V = —X -= X -\- z{ijj\ where uj = x - 



24 i-^ 12 ^ " i+V2^ 

and reduces (30) to the simple linear ODE iozi^^^j + 2zi^ = whose general solution 
z = Cl + C2C0~^ substituted to the ansatz gives the exact solution 
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of equation (30). Applying the inverse hodograph transformation and canceling 
the constant ci by translations with respect to t, we construct the non-Lie solution 

v = v2x — — — 7i (31) 

of equation (29). The solution (31) with C2 = is a Lie solution invariant with re- 
spect to the dilatation operator D = 4tdt + 2xdx + from the maximal Lie 
invariancc algebra of equation (29). However, it is much harder to find this 
solution by the reduction with respect to the operator D. The corresponding 
ansatz v = ^/xz{^o), where u = t~^x^, has a simple form but the reduced ODE 
Alo^z^^ + lo{uj + 4)z^ + 22; — = is nonlinear and complicated. 

This example justifies the observation made by W. Fushchych [7] that "ansatzes 
generated by conditional symmetry operators often reduce an initial nonlinear 
equation to a linear one. As a rule, a Lie reduction does not change the non- 
linear structure of an equation." We can also formulate the more general similar 
observation that a complicated non-Lie ansatz may lead to a simple reduced equa- 
tion while a simple Lie ansatz may give a complicated reduced equation which is 
difficult to be integrated. 

One of the preimages of equation (29) with respect to transformation (4) is the 
equation 

xsm {V2ln x)ut = ( xsiyi'^{V2lnx)ux^ —x ^sin^ 
having the non-Lie exact solution 

[2 r- 3x4 + 24tx2 + C2 

u = \\ — \ cosec V21nx — -. : — 7. . 

V x ' ^ ^' a;4 + 24tx2-c2 

Example 6. Consider the equation from the imaged class (5) 



vt = Vxx-x 2?;"' + — ^ (32) 
16 x'^ 



for the values x > (case 3 of table 1 with m = 3, fc = — |,a2 = ^ and b = — 1). 
It admits the reduction operator of form (27) 

Qj^ = dt + ^ {y2x~iv - x"'^^ \ {^x'^v^ - 3 V2x~3^;2 + x''^V^ dy. 

Usage of the same technique as in the previous example gives the non-Lie exact 
solution of (32) 

1 ^ 1 3t + x2 , , 
i, = -5V2x4— ; . 33 

2 V^15t + x2 +C2 ^ ' 
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Applying the transformation v = ■\/x{hiX'^ + '^)u\,o solution (33), we obtain 
a non-Lie solution of the equation 



4)V 



(^x{blXi+b2X 4)^1(3,^ — ^/x{blX'^ + b2X i ] 

from class (3), where 61 and 62 are arbitrary constants, 6f + 63 7^ 0- 

The equivalence of cases 3 and 4 from table 1 with respect to the point trans- 
formation (7) allows us to use solutions (31) and (33) for finding non-Lie exact 
solutions of the equations 

vt = V,, - x-^P'^'v^ + (-p\' + 1^-' + 2p) V, (34) 
vt = v^^ - x-ief^'-u^ + {-p^x"^ + + l)v (35) 

(case 4 of table 1, where m = 3, (fc, 02) G { (-3, |) , (-|, ^) } and 5 = —1). The 
obtained solutions of (34) and (35) are respectively 

_p 2 , — 3pa;^ - + C2e^^'* 

V = e 2^ V2x — ^ 

px^ — — C2e^'P^ 

and 

Apx"^ - 3 



v = 5\/2e 2^ a;4 



2^(4px2 - 15) + C2e5p* ■ 

7 Appendix 

After reducing an equation from class (5) with the coefficients given by (12), (13) 
(m = 3, fci = 0) by means of operator (14), we need to integrate an ODE of the 
form = Pz'^ + Qz"^ + i?, where P, Q and R arc real constants (see note 1). 
By scale transformations, this equation can be transformed to one from those 
with righthand sides adduced in the fourth column of table 3. The corresponding 
solutions are Jacobian elliptic functions [1,24]. Below 

cn(u;,//) 1 A i \ 

dn(a;,/x) sn(a;,/x) 

nc(w,/x) = — — -, nd{u},n) = —-^ — r-, sc{uj,ii)-- 

sd{uj,n} = — -, cs{uj,n) = — -, ds(a;,/x)- 

dn[LO,n) sc(uj,fi) sd[LO,n) 

The parameter /x is a real number. Without loss of generality ji is supposed to be 
in the closed interval [0, 1] since elliptic functions whose parameter is real can be 
made to depend on elliptic functions whose parameter lies between and 1 [1, §16]. 
If /i is equal to or 1, the Jacobian elliptic functions degenerate to elementary 
ones. 



1 


cd(a; 


,/")' 


sn(a;, 


M) 


cn(a;. 




1 
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Table 3 ( [23]). Relations between values of {P,Q,R) and corresponding 
solutions z{lo) of the ODE = Pz"^ + Qz^ + R. 



p 


Q 




Pz"^ + Qz^ + R 


z{uj) 


A* 




1 
1 




CU. 1 U/ , [Jj J 




2/;2 - 1 








1 


Z // 




( ] -2'| f -2 1 2 -|\ 




1 

J. 






/-i _ ^2\ ( 2 _ 2\ 
[L Z ) Z ) 


TlQI/.l 1(1 

r\n{ ! 1) II] 

\A\^\LU , Ll J 


1 - 


2/^2 - 1 










2-^2 


-1 


(l-z^) ((1-^V-l) 


nd(u;, fi) 


1-/^2 


2-/^2 


1 


(1 + ^2) ((1-^2^^2 + 1) 


sc{lo, fl) 


/.^(/X^-l) 


2/i2_i 


1 


(1 + ((/x2 _ 1)^2 + 1) 


sd{uj, ji) 


1 


2-/^2 




(1 + ^2) (^2 + i_^2) 


cs(a;, /x) 


1 


2^^2-1 


- 1) 


(z2 + ^2) (^2+^2_i) 


ds(a;, /x) 
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